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Introduction 



Soon after the discovery of Quantum Chromodynamics (QCD), it was conjectured that 
at high temperature T the color charge is screened [1] and the corresponding phase of 
matter was named Quark-Gluon Plasma (QGP). It is a special kind of plasma in which 
the electric charges are replaced by the color charges of quarks and gluons, mediating 
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the strong interaction among them. Such a state of matter is expected to exist at extreme 
temperatures, above 150 MeV, or densities, above about 10 times normal nuclear matter 
density. These conditions could existed in the early universe for the first few microseconds 
after the Big Bang or in the interior of neutron stars. 

The aim of the ongoing relativistic heavy-ion collision experiments is to explore the pos- 
sible QGP phase of QCD. The essential difference in heavy-ion collisions over the nucleon- 
nucleon collisions is the dominance of the partonic-level description for essentially all mo- 
mentum scales and over nuclear size distances. In order to describe the produced system 
as a state of matter it is necessary to establish that these nonhadronic degrees of freedom 
form a statistical ensemble. Therefore, the concepts such as temperature, chemical poten- 
tial and flow velocity should apply and the system can be explained by an experimentally 
determined equation of state. In addition experiments should eventually be able to deter- 
mine the physical characteristics of the phase transition, viz., the critical temperature, the 
order of the phase transition, and the speed of sound along with the nature of underlying 
quasi-particles. To date such description is provided by the Lattice QCD (LQCD) calcu- 
lations. Ultimately, one would expect to validate this by characterizing the QGP in terms 
of its experimentally observed properties. The commissioning of Relativistic Heavy Ion 
Collider (RHIC) at BNL and various experiments performed therein have ushered in a new 
era. These experiments have acquired data on Au+Au collisions at various energies, an 
essential p+p baseline data set, and a critical d+Au comparison. The analyses of these 
various system have yielded a rich abundance of results. 

There were non-overlapping talks and overlapping talks during WHEPR The non- over- 
lapping talks on this subject were focussed on the recent developments in LQCD, present 
status of experimental observations and the corresponding theoretical efforts, and certain 
aspects of future experiments, in understanding the properties of QGP produced in heavy- 
ion collisions. Details are reported in these proceedings. 

1. Status of Lattice QCD: Rajiv V. Gavai. 

2. Results from STAR Experiment at RHIC: Bedanga Mohanty. 

3. PHENIX Overviews - Status of QGP: Terry Awes. 

4. Jets in ALICE: J. Guillermo Contreras. 

5. Theoretical Status of QGP: Jean-Y. Ollitrault. 

6. High Energy Photons from Relativistic Heavy Ion Collider: Dinesh K. Srivastava. 

7. On the Structure and Appearance of Quarks Stars: Prashanth Jaikumar. 

Recent numerical LQCD calculations has given us wealth of information on various 
thermodynamic properties at finite temperature and chemical potential. It has now been 
established that there is only a crossover [2] of normal hadronic matter to a state of decon- 
fined quarks and gluons at temperature T c <~ 200 MeV. Moreover, the equation of state [2], 
various susceptibilities [3,4] and transport coefficients [5] have been obtained. It is also 
found that charmonium states remain bound at least up to T <~ 2Tq [6] and the behavior 
of temporal correlators in pseudoscalar and vector channels deviates significantly from the 
free behavior at T ~ 3Tc [7]. These analyses suggested that pseudoscalars and vectors 
resonances may exist above Tc- 
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Robust results from Au+Au at the BNL RHIC experiments have shown collective effects 
known as radial [8] and elliptic [9,10] flows, and a suppression of high-pr hadron spec- 
tra [10,1 1], which could possibly indicate the quenching of light quark and gluon jets [12]. 
The hydrodynamical description of the observed collective flow indicates that the matter 
produced at RHIC behaves like a near-perfect fluid [13]. On the other hand the amount of 
jet quenching might depend on the state of matter of the fireball, i.e., QGP or a hot hadron 
gas. There are extensive theoretical efforts in understanding the effect of the medium on 
jet quenching [14-18] using collisional as well as radiative energy loss since the high en- 
ergy partons traveling through a medium will lose energy owing to the interactions in the 
medium. The measurement in RHIC [19] indicates excess direct photons over the next to 
leading order pQCD processes in RHIC has also been reported and extensive theoretical 
efforts were made to understand these excess photons [20]. In astrophysical part the fo- 
cus was to study and understand the structure of Standard Model at low energy and the 
evolution of the early Universe. 

The working group activity was structured around review talks followed by discussions. 
The emphasis in this working group (IV) was to establish connection between the subject 
matter of the review talks in plenary sessions and some informal talks during discussion 
sessions, which led to identification of relevant physics problems along with efforts to 
work them out. This working group IV was particularly interested in following topics and 
relevant problems: 1 ) Charmonium suppression, 2) Jet quenching and jet identification, 3) 
Hydrodynamic description of elliptic flow, 4) Susceptibilities and speed of sound in QGP, 
5) Possibility of bound states in QGP, 6) ^-production at RHIC, 7)AdS/CFT- QCD/QGP 
and 8) Neutrino emission from crystalline color superconducting quark matter. 

In the next few sections we will briefly report the discussion held on these topics, and 
also the progress made on the relevant physics problems undertaken by the working group 
members. 



1. Charmonium suppression 

L. Ramello 

The latest results from CERN SPS fixed target experiments NA50 and NA60 about heavy 
quarkonia (mostly charmonium) production in p-nucleus collisions at 400 and 450 GeV 
and in In-In, Pb-Pb collisions at 158 GeV/nucleon are presented. For a historical overview 
of charmonium studies at SPS see [21] 

Experiments NA50 and NA60 share the same muon spectrometer which allows to trig- 
ger on dimuons emitted from charmonium and bottomonium decays, as well from the 
Drell-Yan process and open charm associated production. NA50 has three independent 
centrality detectors: an electromagnetic calorimeter, a multiplicity detector and a zero- 
degree calorimeter (ZDC). NA60 retains the ZDC for centrality measurement and has a 
vertex magnet and a silicon pixel vertex spectrometer which allows matching between the 
muon tracks in the muon spectrometer and in the vertex spectrometer. 

The recent NA50 p-A results at 450 and 400 GeV [22] allow (together with NA38 and 
NA3 data at 200 GeV) to determine the J/ip nuclear absorption only from proton induced 
reactions, without reference to S-U results. A value of a absorption^) ^ 'VO = 4.18 ± 0.35 mb 
is obtained. This in turn allows to calculate the expected J/tjj yield in Pb-Pb collisions at 
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158 GeV/nucleon and in S-U collisions at 200 GeV/nucleon. The centrality dependence 
in Pb-Pb is calculated with a Glauber model using the variable L (average path in nuclear 
matter) as a common parameter in p-A, S-U, In-In and Pb-Pb collisions. 

The NA50 (J/ip I Drell-Yan) and (tp' I Drell-Yan) cross-section ratios as a function of 
centrality are obtained from Pb-Pb data samples collected at 158 GeV per nucleon in the 
1998 run and, under improved experimental conditions, in 2000 [23]. 

The J/tp expected production extrapolated from p-A data is then compared to the NA50 
results in Pb-Pb collisions, as well as to the NA38 results for S-U reactions. A departure 
from normal nuclear absorption is observed for mid-central Pb-Pb collisions, with the sup- 
pression increasing with centrality. The three centrality variables give a consistent picture. 
On the other hand, peripheral Pb-Pb data and all S-U data are compatible with normal 
nuclear absorption. The ip' has a significant absorption in p-A collision, which is further 
increased by about a factor 3 both in S-U and Pb-Pb collisions. The J/ip suppression has 
been studied as a function of p^: the anomalous suppression is concentrated at low p^. 
For pt > 3.5 GeV/c, the central to peripheral ratio Rcp becomes almost independent of 
the centrality range. The average square pt (or equivalently the effective temperature) of 
the J/ip first increases with centrality and then saturates in central Pb-Pb collisions. 

The NA60 experiment has studied the centrality dependence of the i/i/j production in 
another system, namely In-In, at 158 GeV per nucleon [24]. The JA/>/Drell-Yan ratio in 
In-In in three centrality bins compares well with NA50 Pb-Pb results, in the sense that 
anomalous suppression is seen also in In-In in the limited range of L accessible to this 
analysis. NA60 has presented more detailed (still preliminary) results by directly compar- 
ing the measured sample with a theoretical distribution expected in the case of pure 
nuclear absorption. The onset of anomalous suppression in In-In is clearly seen in the 
range 80 < N par t < 1 10 (N par t being the number of participating nucleons), with satura- 
tion for larger values of N par t, up to about 200. Comparing the NA60 direct J/ijj analysis 
with NA38/NA50 results, it is evident that the three data sets do not overlap when plotted 
against the L variable, while a rather good overlap is seen when they are plotted against 
the N parf variable, as seen in Fig. 1. 
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Figure 1. The S-U, In-In and Pb-Pb data points do not overlap in the L variable, as 
seen in the left panel. The J/t[> suppression patterns for different interacting systems are 
in fair agreement with the N pa rt variable. 
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The SPS data presented above and the PHENIX results [25] obtained with Au+Au and 
Cu+Cu collisions at RHIC (cm. energy of 200 GeV/nucleon) have been compared with a 
few models [26-29], which attempt to describe data under different assumptions. Models 
have in general been tuned on p-A, S-U and Pb-Pb SPS data and give predictions for In-In 
data and for RHIC Au-Au (Cu-Cu) data. 

Model [26], based on \c suppression in a percolation scenario, predicts an onset of the 
suppression at SPS in Pb-Pb collisions at N part ~ 125, which agrees with NA50 data, and 
an onset in In-In collisions at N port ~ 140, which disagrees with NA60 data (the onset in 
data is at N part <~ 90). 

Model [27] contemplates charmonium suppression and regeneration in both QGP and 
Hadron Gas phases. At SPS only a small amount of regeneration is needed, and both 
NA50 and NA60 data are well described, except for In-In data at N part > 160 which are 
underpredicted. At RHIC 200 GeV cm. energy in Au-Au collision a substantial amount 
of regeneration is needed to describe data. 

Model [28] includes nuclear absorption (a a bs = 4.5 mb) and comover absorption (er co = 
0.65 mb). SPS data for Pb-Pb collisions are well described, while those for In-In collision 
are slightly underpredicted: the ratio Measured/Expected (expected from normal nuclear 
absorption) for semi-central and central collisions is 0.80 in data with no centrality depen- 
dence, 0.70 and decreasing in the model. The R^ ratio at RHIC for Au-Au and Cu-Cu 
collisions at N part ~ 300 is largely underpredicted (0.10 vs. 0.35). 

Model [29] has in fact two scenarios: the QGP scenario underpredicts the Au-Au 
PHENIX data, while the Statistical Coalescence Model (valid for N part > 100) predicts a 
flat behaviour vs. N part and roughly agrees with PHENIX dN(J/?/>)/dy. 

2. Jets 

2.1 Jets in ALICE 
J.G.Contreras 

Jets, at the parton level, are cascades of quarks and gluons produced by a fast moving 
parton in a process called fragmentation. In the ALICE experiment, the fast partons are 
created as a product of a hard scattering, where the incoming partons are provided by the 
collision of two extremely relativistic protons or nuclei. At the detector level, jets are the 
highly collimated bundle of particles produced by the hadronization process of the partonic 
jet. 

The fragmentation properties of a parton jet are influenced by its interaction with the 
medium, where the parton loses its energy in the QGP. This phenomenom is called jet 
quenching and it has been advocated to be an ideal probe of the QGP state created by 
the collision of two heavy ions at asymptotic energies [12]. Partonic jets from a hard 
scattering are created before the QGP in the collision of two heavy ions. Hard scatterings 
and the evolution of fragmentation in vacuum can be quite reliably computed within pQCD. 
Any changes in the fragmentation process in the presence of QGP can be determined, in 
principle, from the measurement of jet quenching. There are some caveats in this approach. 
The experiments measure jets after hadronization, the effect of which can not be calculated 
a priori, inhibiting a unique determination of the effect of QGP on the partons. However, 
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there are models which assume (with certain confidence) that the influence of hadronization 
is small in certain regions of phase space. Also, lacking detailed understanding of the 
QGP state, the understanding of the interaction of a parton jet with a QGP is still in an 
evolutionary phase. 

While there have been attempts to understand the effect of the medium on the jet quench- 
ing due to the radiative energy loss [14-17], and recently also due to the collisional energy 
loss [18]. Furthermore, only the first approximation to radiative energy loss has been used, 
but a lot of theoretical work is under progress in this area. Not withstanding these limi- 
tations, observation of jet quenching in RHIC data [11] has been used to 'announce' the 
creation of a QGP. 

The experiments at RHIC have shown that there is jet quenching and that it is definitely 
a final state effect. It must be mentioned here that there is an inherent bias in leading 
particle analysis, which can be best eliminated by studying the complete jets. Further, it is 
desirable to study the QCD evolution of jet quenching. This can be best done by studying 
evolution of jet quenching with the hardness of the scale by comparing jet shapes in a wide 
kinematic domain. These steps lead to the ALICE experiment at the Large Hadron Collider 
at CERN; to new problems, and also to new opportunities. 

The physics reach of studying jets in the ALICE experiment has been reported in [30] 
and a brief overview is discussed here. 

It is expected that under nominal conditions of LHC running, several tens of jets with 
energies below 15 GeV will be produced in each Pb-Pb event within the acceptance of 
ALICE. It is also expected that there will be jet(s) of 250 GeV for one month of data taking. 
The large charge particle multiplicity envisaged in the heavy ion collision, in contrast to 
pp collision, will lead to several hundreds of GeV of energy deposition in one unit in the 
azimuth-rapidity space. The fluctuations in the average energy in a given cone are of the 
order of several tens of GeV and form a formidable background to observe jets. 

The strategy planned for the ALICE experiment is to use smaller cones, and to search 
for the jets using particles whose energy lies above a given threshold and to iteratively 
subtract the background using a modified jet cone algorithm [31]. It has been shown that 
above 50-60 GeV of energy, jets and its spectra can be reliably reconstructed using this 
modified jet algorithm on the charged particles produced in the interaction and measured 
by the TPC. 

Jet quenching affects certain properties of the jet like amount of radiation outside a given 
cone, jet heating and the fragmentation function [15]. The various variables have relative 
merit e.g. one variable (jet heating) does not evolve with the hardness of the interaction 
and can not be used to test the QCD evolution of jet quenching, while remaining a good 
test of jet quenching at all hardness scales. In contrast, fragmentation function of jets offers 
a window to study the evolution of jet quenching with the hardness of the interaction. 

In summary, jet shapes are ideal testbeds to study jet quenching phenomena. To use 
them it is mandatory to measure very low energy particles and to understand with great 
precision the contribution from the underlying event and its fluctuation. ALICE fulfills 
all the requirements to perform these demanding measurements and extract from them the 
properties of jet quenching and thus increase our understanding of the QGP. 
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2.2 Flow Coefficients and Jet Characteristics in Heavy Ion Collisions: New Methods of Jet 
Identification 



S.C.Phatak 



Identification of jets in heavy ion collisions is an important and challenging problem be- 
cause modification of jet properties is expected to give information on possible formation 
of quark-gluon plasma during the collision process [14-18]. Jet quenching has already 
been observed [11], although the observation is somewhat indirect. In particular, these 
studies are basically the correlation studies between the energetic hadrons which are ex- 
pected to be the hadrons produced in a jet. Detecting and identifying jets in heavy ion 
collisions is an essential prerequisite to study jet quenching. 

A new method for jet identification has been developed [32]. The method is based on 
the fact that the flow or Fourier coefficients for events containing jets have typical structure 
enabling identification of jet events and determining its opening angle and the number of 
particles. The even coefficients in back-to-back jets are observed to be larger. 

The flow coefficients are Fourier coefficients of the azimuthal distribution of particles 
produced in heavy ion collisions and are given by 

i iPr = J d(t>ld(t>2dpTldpT2PTlPT2P{<t>l)P{<t>2) cosm(<f>i - &) (1) 

Given the experimental distribution of particles, the flow coefficients are given by pT = 
~W Si j Pt.iPtj cos((f>i — <pj) where N is the number of particles in the event and px,i 
is the transverse momentum of i th particle. For uniformly distributed particles, all flow 
coefficients ( except for m=0 ) vanish where as for ^-function distribution all coefficients 
are unity. Thus, for an event with a jet having a number of particles in a small azimuthal 
angle and the rest of the particles uniformly distributed, v m 's are expected to be abnormally 
large. This can be used to identify a jet and to determine the jet properties. We have shown 
that [33], if an event has jet particles distributed in a narrow cone of azimuthal angle, we 

have^ = Nj< ^i > 1 — m 2 a 2 +0(m A ) where a is the variance of the (^-distribution 
and < pt > is the average transverse momentum carried by a hadron in the jet. Thus 
Nj < pt > is the transverse momentum of the jet. 

These expressions clearly suggest a method of obtaining jet properties from the flow 
coefficient. A linear fit to v^ VT , plotted as a function of m 2 would yield the number of jet 
particles and jet px from the intercept on y axis and a from the slope. 

Let us now consider the case of two back-to-back jets in a background of uniformly 
distributed particles. This case is of interest because we expect that a hard parton scattering 
would produce such jets having equal and opposite jet momenta. We expect that quenching 
of one of the jets would broaden one of the jets and/or produce more jet particles. Thus, 
the characteristics of the two back-to-backjets would be different. Further, in an extreme 
situation, the fast moving parton of one of the jets may be completely absorbed in the 
medium leading to removal of one of the jets. Assuming that the jet particles are distributed 
uniformly in the jet cone of the respective jets, we get 



2 J_ 

V m,p T jy2 



< PT,1 > JO(^T^) + (-I)™ < PT.2 > JO(^T^) 



(2) 
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Here Nji are the jet particles in the ith jet, A<j>i are the corresponding opening angles and 
< Pt.i > the average transverse momenta of the respective jet particles. The presence of 
(— l) m factor in the second term above implies odd-even staggering in v^'s. In particular, 
odd coefficients vanish for identical jets. This property would be useful in estimating the 
possible quenching of one of the jets. The method also works when transverse momentum 
is not measured and only the azimuthal angle of the particles is known. 

The method was tested using simulations. Jet particles were simulated by adding certain 
number of particles distributed uniformly in a jet cone of A(j> and having exponentially 
decreasing px distribution and added to HIJING events [34] adapted for LHC energies. 
Different values of number of jet particles and jet opening angles were used for analysis. 
In the absence of any jet particles, the method yields d^ Pt values smaller than 0.05 Gev 2 . 
In case of the data with jet only, the method outlined above is able to extract correctly the 
number of jet particles, jet transverse momentum and jet opening angle. 

3. Dissipative relativistic hydrodynamics 

R.S. Bhalerao 

Certain global features of the RHIC data on ultrarelativistic heavy-ion collisions have 
been successfully explained in the framework of the boost-invariant ideal hydrodynamics: 
The main characteristics of the observed elliptic flow iv 2 ) are reasonably well described by 
ideal fluid dynamics, while requiring unreasonably large cross sections in transport models. 
The ability of the former to reproduce both the elliptic flow and single-particle spectra for 
measured hadrons with pr < 2 GeV/c near midrapidity in minimum-bias collisions is 
considered a significant finding at RHIC. The dependence of the flow pattern on hadron 
masses further supports the hydrodynamical picture. 

However, the short time scale for equilibration is difficult to account for microscopi- 
cally [35]. (It has been argued recently that plasma instabilities may provide a mechanism 
for fast thermalization.) Moreover, several features of the data clearly signal the break- 
down of the ideal hydrodynamical description [36]: Hydrodynamic models seem to work 
for minimum-bias data but not for centrality-selected pion and antiproton data. Indeed, 
the pion data presented in Fig. 36 of [37] show that the measured v 2 in central collisions 
overshoots the value obtained by a hydrodynamical computation, w 2 (data) > v 2 (hydro), 
contrary to the expectation that the v 2 measurements saturate the hydrodynamical limit. 
See also Fig. 6 of [37] where the agreement between v 2 (data) and v 2 (hydro) is far from 
satisfactory. This indicates incomplete equilibration and/or importance of dissipative ef- 
fects. 

Hydrodynamic calculations have shown that v 2 is independent of the system size for a 
given shape of the collision zone. This is a consequence of scale invariance of ideal hy- 
drodynamics. However, incomplete equilibration breaks this scale invariance making v 2 
depend upon the system size. Certain dimensionless numbers like Mach number, Knud- 
sen number and Reynolds number can be used to characterize the motion of the fluid, in 
particular its compressibility, degree of thermalization and importance of viscous effects. 
Estimating these numbers for Au-Au collisions at RHIC calls in question the assumption of 
local equilibrium and applicability of ideal hydrodynamics, and at the same time points to 
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the need for dissipative relativistic hydrodynamic description of nucleus-nucleus collisions 
atRHIC. 



4. Susceptibilities and speed of sound in QGP in PNJL model 

Sanjay K. Ghosh, Tamal K. Mukherjee, Munshi G. Mustafa and Rajarshi Ray 

Susceptibility is the response of the system to an externally applied force. Quark number 
susceptibilities are the response of the quark number density with the variation of chemical 
potential. Lattice data on QCD thermodynamics, particularly recent study of higher order 
susceptibilities [3,4] have provided valuable information and new insight about the proper- 
ties of matter produced around and above the critical temperature, Tc- Perturbative QCD 
calculations fail to describe the details of these results. There are some model calculations 
to understand the physical picture of the structures in higher order susceptibilities produced 
by those LQCD data. The hadron gas resonance model [38] describes the data well below 
Tc but fails for T > Tc- The recently proposed [39] scenario of colored bound states also 
compares the data. However, their comparison is still not completely satisfactory. In this 
work 1 we try to understand the features and structures of the susceptibilities produced by 
LQCD within the model described below. 

The thermal average of the Polyakov-loop can be considered to be the order parameter 
for deconfinement transition [40]. Hence a judicious use of the Polyakov-loop in effective 
models may prove to be of great advantage. On the other hand, there are QCD inspired 
models, predominantly, the NJL model, has given rise to an interesting phase diagram in 
temperature and chemical potentials. 

We study [4 1 ] some of the thermodynamic properties of strongly interacting matter using 
the Polyakov-loop + Nambu-Jona-Lasinio (PNJL) model [42]. The motivation behind the 
PNJL model is to couple the chiral and deconfinement order parameters inside a single 
framework. We have computed the EOS, the quark number susceptibilities, the specific 
heat Cy, the speed of sound (basically its square, v 2 s ), and the conformal measure C = 
A /e, where A = e — 3P is the interaction measure and e and P are respectively the energy 
density and pressure of strongly interacting matter. Comparisons of these quantities with 
those obtained on the lattice were made. 

We use the parameterization of the PNJL model as given in Ref. [43]. For our results we 
shall present the quantities as a function of temperature in units of a crossover temperature 
which is taken to be T c = 227 McV. 

We found the pressure to grow from almost zero at low temperatures to about 90% of its 
ideal gas value at 2.5T C . Left panel in Fig. 2 shows the variation of the quark number sus- 
ceptibility (QNS), which is the second order coefficient (c 2 ) of reduced pressure (P/T 4 ), 
with T/T c . This shows an order parameterlike behavior. At higher temperatures the c 2 



1 The authors are thankful to Rajiv Gavai for useful discussion during WHEPP-09 where this work 
was initiated. Full details of this work was reported in hep-ph/0603050 and published in Phys. Rev. 
D73, 114007 (2006). 
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Figure 2. The QNS (left) and the forth order coefficient of reduced pressure (right) as 
a function of T/T c . Symbols are Lattice data. Arrows on the right indicate ideal gas 
values. 
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T/Tc , T/Tc 

Figure 3. CV/T 3 and 4e/T 4 (left) and Squared velocity of sound vl and conformal 
measure C = A/e (right) as function of T/T c . The arrow on the right shows the ideal 
gas value. 



reaches almost 85 % of its ideal gas value, consistent with Lattice data. Similar behavior 
has been observed in another model study [44] using density dependent quark mass model. 

The fourth order derivative C4, which can then be thought of as the "susceptibility" of 
C2 shows a peak at T = T c (right panel of Fig. 2). Near the transition temperature T c , the 
effective model should work well and we observe that the structure of C4 is quite consistent 
with Lattice data. Just above T c however, there is a significant difference between our 
results of C4 and that of Ref. [3]. Note that in the SB limit both C2 and C4 have only 
fermionic contributions. We expect that because the coupling strength is still large in this 
temperature regime it is unlikely that C4 should go to the SB limit within T < 2.5T C . 
Moreover, the quark masses used in Ref. [3] is considerably large (m/T — 0.4) to expect 
fermionic observables to go to the SB limit. However, it is possible that our overestimation 
is due to the use of mean field approximation. 

Left panel of Fig. 3, shows that Cy grows with increasing temperature and reaches a peak 
at T c . For comparison, we have also plotted the values of 4e/T 4 , at which the specific heat 
is expected to coincide for a conformal gas. 

We plot the speed of sound and conformal measure in right panel of Fig. 3. The value 
of P/e matches with that of v 2 for T < T c and also goes close again above 2.5T C . But in 
between these two limits the v 2 is distinctly greater than PJ e. Thus, C would go to zero 



10 



much faster if we replace P j e by v 2 for computing C. Near T c the v 2 has minimum which 
goes close to 0.08. 

The values for the conformal measure C also closely resembles the Lattice data of Ref. 
[45]. The dip at temperatures less than T c is prominent in both the cases. At even lower 
temperatures we find C to increase. For a nonrelativistic ideal gas, the ratio of P/e should 
go to zero, and thus C should then go to 1 . On the other hand at high temperatures either 
an ideal gas or a conformal behavior should be recovered for which C should go to zero. 

5. Possibility of bound states in QGP 

Sanjay K. Ghosh, Munshi G. Mustafa and Rajarshi Ray 

Recent LQCD calculations [6,7] suggest that pseudoscalars and vectors resonances may 
exist above Tc- Color bound states of parton at rest have been claimed by analyzing LQCD 
data [46]. The situation, however, changes if the test charge, Q a , is in motion relative to the 
heat bath. The motion of the particle fixes the direction in space and spherical symmetry 
of the problem reduces to axial symmetry. This implies the loss of spherical symmetry 
of the Debye screening cloud around the moving test charge resulting in a wake in the 
induced charge due to dynamical screening [47]. The negative minimum found in the 
wake potential indicates a induced space charge density of opposite sign. Thus, a particle 
moving relative to a particle in the induced space charge density would constitute a dipole 
oriented along the direction of motion. We briefly describe the potential due to such dipole 
interaction in QGR details of which can be seen in Ref. [47]. 

The spatial distribution of the scaled 2 potential in cylindrical coordinates are displayed 
in Fig. 4 for two velocities. The left panel of Fig. 4 corresponds to v = 0.55c where 
as the right panel is for v = 0.99c. The dipole potential shows the usual singularity 
of the screening potential at r = (z — and p = 0) and a completely symmetric 
behavior along with a pronounced negative minimum in the (p — z) plane. This is due to 
the compensating effects between the electric and the magnetic interactions. For details 
of electric and magnetic interactions, readers are referred to [47]. With the increase of v 
the relative strength of the dipole potential grows strongly. For v — 0.99c a substantial 
repulsive interaction has grown in the transverse plane, i.e., in the direction of p but at 
z = 0, becoming responsible for a vertical split of the minimum in the (p - z) plane. 
However, the general form of the potential mostly resembles the Lennard-Jones type with 
a pronounced repulsive as well attractive part. 

In QCD the interaction between color charges in various channels is either attractive or 
repulsive. A quark and an antiquark yield the sum of irreducible color representations 3 : 
3 (g) 3 = 1 © 8, with the interaction strength of the color singlet representation is -16/3 
(attractive) whereas that of the color octet channel is 2/3 (repulsive). Similarly, a two 



2 The potential is scaled with the screening mass mo as well as with the interaction strength, Q a Q b . 
So, the details of the potential will depend on the temperature, the strong coupling constant and the 
sign of the interaction strengthi [47]. 

3 Their strength of the color interaction [46] can be calculated using SU (3) color group 
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Figure 4. Left panel: Spatial distribution of the scaled dipole potential with respect 
Debye screening mass mo for v = 0.55c. Right panel: Same as left panel for 
v = 0.99c. 

quark state corresponds to the sum of the irreducible color representations: 3 ® 3 = 3 © 6, 
where the antisymmetric color triplet is attractive (—8/3) giving rise to possible bound 
states [46]. The symmetric color sextet channel, on the other hand, is repulsive (4/3). 
Color bound states (i.e. diquarks) of partons at rest have been claimed by analyzing LQCD 
data [46]. The situation is different when partons are in motion. The dipole potentials 
along the direction of propagation and also normal to it have both attractive and repulsive 
parts, similar to the Lennard-Jones form. So, all the attractive channels or the repulsive 
channels in the static case get inverted due to the two comoving partons constituting a 
dipole in the QGP. This could lead to dissociation of bound states (or to resonance states) 
as well formation of color bound states in the QGP. 

Within our model such bound states as well other colored binary states in the QGP can 
experience different potentials along the dipole direction and the direction normal to it. 
Along the direction of motion binary states which were bound in QGP may become reso- 
nance states or dissociate beyond Tq- Similarly, those colored states which were not bound 
initially in the QGP, may transform into bound states. There are some long distance corre- 
lations among partons in the QGP, which could indicate the appearance/disappearance of 
binary states in the QGP beyond Tq ■ The temperature up to which they survive need fur- 
ther analysis of the bound states properties in detail, which was initiated during WHEPP-09 
and some progress has been made. The complete work will be reported elsewhere. 

6. ^-production at RHIC 

Jajati K. Nayak, Jan-e Alam, Bedangadas Mohanty, Pradip Roy, Abhee K. Dutt-Mazumder 

Among the promising signals to analyse the matter formed in Relativistic Heavy Ion Col- 
lision, the study of strangeness is important one. The production and evolution of strange 
particles give good information about the formation of QGP. <f> meson is a bound state of 
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strange quark(s) and anti strange quark (s). The interaction of <f> meson with nuclear matter 
is suppressed according to Okubo-Zweig- Izuka(OZI) rule. <fr meson after its production 
during hadronisation suffers less rescattering with hadronic matter. So it gives thermody- 
namic information of the state of matter during hadronisation stage. From the <j> spectra at 
^=200 GeV (RHIC energy) we extract the information about critical temperature T c of 
the quark-hadron phase transition. 
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Figure 5. Evolution of strange particles and temperature when T c =170 MeV. rii/r\ q 
is the ratio of no. density of i-type particle. 

We assume a first order quark-hadron phase transition for the RHIC energy -^=200 
GeV. In QGP phase s and s quarks are produced mainly by gluon gluon fusion and an- 
nihilation of light quarks and anti-quarks. These s and s forms cf> through hadronisation 
process. This production is not OZI suppressed. We expect excess cf> mesons if QGP is 
formed in the initial state of Heavy Ion Collisions as compared to a hadronic initial state. 
The production of <fi from sources other than plasma is very small. <j> spectra is experi- 
mentally reconstructed from Kaons since it decays to K + and K~ . Although K + and 
K ~ carry some scattering effect still <j> meson spectra gives reliable information of the 
thermodynamic state of matter during its formation. 

We calculate the ratio of the multiplicity of <f> meson to the total multiplicity at mid 
rapidity of RHIC experiment from thermal model. The ratio depends on the critical tem- 
perature T c and the effective degeneracy g e ff- We took these value from lattice QCD 
calculation. We found the measured value to be more than the calculated one. This is ex- 
pressed as a over production or enhancement factor 70. To know where and when this over 
production occurs we adopted Boltzmann equation and solved for the evolution of strange 
particles(e.g, s,s, K + ,K~ etc) in QGP, mixed and hadron phases. With proper initial con- 
ditions from hadron spectra we solved the Boltzmann equation and temperature evolution 
equation and found the over production factor 70 r| f+ =1.63 (rjc+ being the ratio of non- 
equilibrium to equilibrium density of K + ) when T c =170 MeV. The production of K + is 
restricted by a parameter 8. Maximum half of the strange quarks can go to the K + . This 
is expressed as (5=0.5. So this puts the limit on the critical temperature i.e, T c > 170 MeV. 
Similarly a limit on g e ff also has been put [48]. 
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7. AdS/CFT ► QCD/QGP ? 

Balram Rai 

The idea that QGP at RHIC, T <~ 2Tc, seems to be in a strongly coupled regime (interac- 
tion energy > > kinetic energy) and has presently been attracted a lot of attention : 

1 . The equation of states obtained in LQCD deviates from ideal gas behavior. Com- 
bining the LQCD data on quasiparticle masses and interparticle potential one indeed 
finds [46] a lot of bound states and resonances above Tc- This could be a explana- 
tion for charmonium states remaining bound up to (2-3)7b, as was directly observed 
on lattice [6]. However, resonances above Tc due to marginal bound states give rise 
to large cross-section and small mean free path leading to collective flow. 

2. Collective phenomena observed at RHIC lead to view QGP as a near perfect 
fluid [13]. The deviation of elliptic flow data from the prediction of hydrodynamic 
calculations only occur at pr ~ 1.5 — 2 GeV/c and this leads to viscosity to entropy 
ratio, rj/s ~ 0. 1-0.2, which is more than order of magnitude less than that in pQCD. 

3. The charm diffusion constant [49] deduced from the single electron data and elliptic 
flow is also an order of magnitude less than pQCD estimates. 

4. The interaction parameter T = (potential energy) /T is obviously not small for 
QGP [50]. At such T the classical strongly coupled plasma is a good liquid. 

As seen the pQCD calculations of transport coefficients are not reliable around and above 
Tc, which is believed to be in strong coupling regime. On the other hand the LQCD is a 
rigorous calculational method applicable in hot, strongly coupled, gauge theory. Because 
it is formulated in Euclidean space, it is well suited for computing static thermodynamic 
quantities and less well-suited for transport coefficients, or dynamical processes of any sort. 
Complementary nonperturbative techniques are thus desirable. One such technique is the 
Anti-de Sitter/Conformal Field Theory (AdS/CFT) correspondence, which maps nonper- 
turbative problems in certain hot strongly coupled gauge theories onto calculable problems 
in a dual gravity theory [51]. This method has been applied to calculate the shear viscos- 
ity [52] in Af — 4 supersymmetric Yang-Mills (SYM) theory which gives an upper bound 
of rj/s <~ 1/47T, close to the RHIC value. Also effort [53] was made to compute certain 
diffusion coefficients in the same spirit. The ratio of pressure to that of Stefan-Boltzmann 
limit in J\f = 4 SYM is remarkably close to the corresponding ratio in LQCD at tempera- 
tures a few times Tc where it is strongly coupled [54]. 

Usually the properties of Af = 4 SYM theory are completely different from QCD. The 
former one is a conformal theory with no particle spectrum or 5-matrix, whereas the latter 
one is confining theory with a realistic particle interpretation. In view of this members 
of this working group were interested to know the basics of AdS/CFT correspondence to 
supersymmetric gauge theories in strong coupling limit and their connection to QCD/QGP. 
Balram Rai delivered two comprehensive lectures on this topic in a very elementary level, 
which are discussed below in brief. 

• Af = 4 SYM is a conformally invariant theory with two parameters: the rank of 
the gauge group Nc and the t'Hooft coupling A = gyM^c with <jym is Yang- 
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Mills coupling. Its on-shell field content includes eight bosonic and eight fermionic 
degrees of freedom, all in color adjoint representation. 

• Basic argument leading to this correspondence: consider in type 1IB string theory, a 
configuration of Nc £>3-branes stacked on top of each other. The appropriate theory 
living on the branes is M = 4 SYM theory. On the other hand if Nc is large, the 
stack of Z?3-branes has large tension, which curves space-time. In the limit of large 
t'Hooft coupling A, the brane geometry has small curvature and can be described by 
supergravity. Therefore, one can have two descriptions of the same physics in terms 
of strongly coupled gauge theory on the branes and classical gravity on a certain 
background. 

• AdS/CFT conjecture (or gauge/string duality) states that this theory is exactly equiv- 
alent to type IIB string theory in an AdS$ x S 5 gravitational background, where 
AdS$ is a five dimensional anti-de Sitter space and S 5 is five dimensional sphere. 
At large N c and large A, the string theory can be approximated by classical type IIB 
supergravity. 

• The above approximation permits nonperturbative calculations in quantum field the- 
ory mapped into problems in classical general relativity. In this context, raising the 
temperature of the gauge theory corresponds to a black hole (or a black brane) into 
the center of AdS$. According to AdS/CFT, the Hawking temperature of the black 
hole becomes the temperature of the gauge theory. 

• Even at T = 0, the M = 4 SYM theory is a conformal one whereas the QCD is 
confining theory, but T / both theories describe hot, nonabelian plasmas with 
Debye screening, finite spatial correlation length, and qualitatively similar hydrody- 
namic behavior. The major differences is that all the excitations (quarks, gluons and 
scalars) in N = 4 SYM plasma are in adjoint representation, while hot QCD plasma 
has quarks in fundamental and gluons in adjoint representation. This leads to an 
assumption that the fundamental representation fields have negligible influence on 
bulk properties of the plasma. One may view the quarks as test particles which serve 
as probes of dynamical processes in the background of J\f = 4 plasma. 

• As found the shear viscosity, the thermodynamic quantities such as pressure, energy 
density, entropy etc. computed in Af SYM in strong coupling limits (large Nc and 
large A) are insensitive to details of the plasma composition or the precise interaction 
strength. These values are remarkably close to those properties of QGP observed in 
RHIC and obtained LQCD. 

• These observations have led to think a connection of AdS/CFT gravity dual to RHIC 
collisions ! 



8. Neutrino emission rates in crystalline color superconducting quark matter 

Prashanth Jaikumar, Hiranmaya Mishra and Andreas Schmitt 

The existence of deconfined quark matter inside compact stars can be verified if the 
stellar cooling rate is found to be substantially different than for neutron stars, and can 



15 



be attributed directly to the neutrino emission rates of (superconducting) quark matter. 
This requires 4 1) identifying the ground state of dense QCD at intermediately large quark 
chemical potentials {p, q ~ 300 — 400 MeV) that characterize the interior of neutron stars, 
and 2) computing the neutrino emission rate and specific heat of this phase to determine its 
cooling behavior. 

One possible manifestation of the diquark phase at intermediate densities, where the 
strange quark mass is large, is a crystalline color superconductor in which quarks with 
different Fermi surfaces pair at non-zero momentum, resulting in an inhomogeneous but 
spatially periodic order parameter [55]. This phase spontaneously breaks translation and 
rotational symmetries, and the free energy of the system is minimized when the gap varies 
spatially in accordance with the residual discrete symmetries of this phase. This is the QCD 
incarnation of the LOFF (Larkin-Ovchinnikov-Fulde-Ferrell) phase in electronic spin sys- 
tems. 

The presence of color superconducting phases is expected to alter the cooling behavior of 
neutron stars. The equation that determines the cooling of an isothermal neutron star of 
volume V in global heat balance is 

dT 

c v — = -L v = -e v V (3) 
dt 

where T is temperature, t is time, e is the neutrino emissivity and c v the specific heat at 
constant volume V of the quark phase. 

Minimization of the thermodynamic potential leads to the nonisotropic gap equation [56] 



2G f , 2 sin 2 | 

—3 / dp 2 (4) 

~ ' ' ^(\q + p\ + \q-p\-fl) 2 +4A 2 sin 2 f 

where the momentum integration runs over the pairing region defined by the domain 
{V : p|-Ei(p) > 7 E 2 (p) > 0}. Here, the angle is the angle between the up quark 
momentum (q + p) and the down quark momentum (q — p). Further, the quasi particle 
energies are given as [57] 

E 1 (p) = +6 f i+^(\p + q\-\p-q\) 



and 



+\J (\Q +p\ + \q-p\- 2m) 2 + 4A 2 sin 2 (i/3) (5) 



E 2 (p) = -Sn- ^(\p + q\ -\p-q\ 



+ 2 Y (l« +P\ + \Q-P\- 2 M) + 4A2 sin ^ 



4 This work was initiated during WHEPP-09 and substantial progress has already been made. 
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In the above, /l = (/i u + /Ltd)/2 and 8(x = \i& — \i u is the average and the difference 
respectively of the chemical potentials of the two condensing quarks. Thus, crystalline 
superconductivity is characterized by dispersion relations which vary with direction of 
momentum, yielding gaps which vary from zero up to a maximum of A. 

For the purpose of neutrino emissivity calculations, it is the nodal surfaces described by 
£i=0 and E2 = that make the dominant contributions to the phase space integral in 
Eqn.(7). Thus one needs to know the dispersion relation near the 'blocking' regions B u , Bd 
in Fig. 6, whose boundary can be specified by solving Ei(p) — for the angle between 
the condensate momentum q and the relative momentum between the pairing quarks p. 





Figure 6. Allowed (unshaded) and blocked (shaded) regions for pairing between up 
and down quarks in the LOFF phase. The 2-dimensional projection of the Fermi spheres 
for the two flavors of quark are shown displaced by an amount 2q. The dispersion 
Ei(E2) vanishes on the boundary of the lower (upper) banana. These surfaces make 
the principal contributions to the specific heat and neutrino emissivity. Particle-hole 
pairing occurs on the boundary of the elliptical region but is clearly disfavored at low 
temperatures. 



The neutrino emissivity is calculated from 



n 



<Pp r 

- = -JE(p,)=v (2tt) 3 
x4>d(pdfBi[> u (p u ) 2 



E v Wn(p d )(l - n(pu))(l - n(p e )) 



where i= 1..4 represents d, u, e, v respectively. 



W = 



(2ir) 2 5 A {p d -p u -p e 



(7) 



(8) 



with |M| 2 = QAG 2 F cos 2 6 c (pd-p v ){pu-Pe) being the matrix element for the weak inter- 
action process. The ip's are Bogoliubov coefficients for the (du) quasiparticles. 
The specific heat can be calculated from 
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, - 1 f d3 p E ^ 2 m 

T 2 J (2tt) 3 (e £ (p)/ T + l)(e- £ W/ T + 1) W 

Since we are interested in the low-temperature behavior of the specific heat, we can 
approximate 

c ^ h I wf E ^ 2e ~ E(p)/T (10) 

The temperature dependence of the specific heat can be classified with respect to the 
behavior of the dispersion relations in the vicinity of its zeros and the dimension d = 
(point), 1 (line) and 2 (surface) of the sub-manifold where the dispersion vanishes. Using 
the above formula, we can deduce the results for the temperature dependence of the specific 
heat in the following table: 



Cy 





1 


2 


linear 




/iT 2 




quadratic 


flT 2 


£3/2 T 3/2 


p6/2 T l/2 



where constants of proportionality of order one have been omitted. These relations in- 
dicate that the circle of nodes at the tips of the two banana-shaped regions contribute the 
most to the specific heat since p, ^> T. The point nodes contribute less to the specific 
heat but will dominate the contribution to the neutrino emissivity since their phase space 
weights are larger. The computation for the emissivity is complicated by the need to per- 
form a bounded phase space integration excluding the blocking regions, and this will be 
attempted numerically. 
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